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An experimentally feasible scheme for generating a 12-qubit hyperentangled state via 
four photons, entangled in polarization, frequency and spatial mode, is proposed. We 
study the nature of quantum non-locality of this hyperentangled state by evaluating 
its violation degree to a Bell-type inequality, and find that the result agrees well with 
quantum mechanics prediction while extremely contradicts to the local realism con- 
straint. 
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1 Introduction 

Qubit is the basic information unit of quantum computer. Hence, the generation of multi- 
qubit entangled states keeps as an essential task in quantum information processing and 
transmission [1]. Quantum entanglement, especially multipartite entanglement, has many 
important applications, such as measurement-based quantum computing [2], secure super- 
dense coding [3], teleportation [4], entanglement purification [5j and quantum cryptography 
[6]. However, the experimental realization of multipartite entanglement is still a significant 
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challenge. People realized that the manipulation of more than six photons without quantum 
storage turns out to be an insurmountable obstacle with present technology [7]. However, 
the hyperentanglement(HE), which in practice remains in many physical systems, tends to 
be an effective and practical way [H [9l [10] to manipulate more qubits effectively, namely 
making the quanta, e.g. photons, to be entangled simultaneously in multiple degrees of 
freedom (DOFs). Furthermore, HE states are much less affected by decoherence and are 
particular important in the realization of even more challenging objectives in comparison 
with the normal entangled states. 

Recently, two-photon four-qubit states [HI [121 [13] and two-photon six-qubit states [T^[T5] 
were experimentally realized, and found to offer significant enhancement of the channel 
capacity in superdense coding [161 [U]) to construct efficiently the fourqubit two-photon 
cluster states [18], and to realize successfully the basic quantum computation algorithms in 
the one-way model [HI [19]. More recently, Jian-Wei Pan et al. generated in experiment 
the hyperentangled ten-qubit photonic Schrodinger cat state by exploiting polarization and 
momentum of five photons [7j . In this article, we propose an experimentally practical scheme 
for generating a four-photon twelve-qubit hyperentangled state entangled in three DOFs, i.e. 
in polarization, frequency and spatial mode. 

Nonlocality is a unique feature of Quantum Mechanics(QM) differing from the local 
realism(LR) theory. People found that the HE states can enhance the violation of LR [20] . 
We find in this work that the twelve-qubit hyperentangled state exhibits extremely large 
violation of the LR constraint by evaluating a Bell-type inequality. Furthermore, we show 
generally that the more the entangled degrees of freedom are, the stronger the violation of 
Bell-type inequality will be 

The structure of the paper goes as follows: in section 2, an experimentally feasible scheme 
for the production and measurement of a four-photon twelve-qubit hyperentangled state 
is proposed. In section 3, we discuss the quantum characteristics of the hyperentangled 
state(mathematic details are presented in Appendix). A brief summary and concluding 
remarks are given in section 4. 
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2 Preparation scheme for hyperentangled state 



The hyperentangled state to be generated in our scheme takes the following form up to a 
normalization constant: 



= {\HHHH) + \VVVV)) (g){\u,U2UJ,U2) + \uj2U^U2U^)) (g){\a[b[c[d[) + ^b'^c'^d'^)) , (1) 



where H and V denote horizontal and vertical polarization, ooi and 002 signify different fre- 
quency and a[, b[, c[, d[, a'2, b'2, c'2, d'2 label eight different spatial modes. The state ([I]) shows 
maximal entanglement among all photons' degrees of freedom in polarization, frequency and 
spatial modes. This state is a Schroinger cat state, and can also be called as Greenberger- 
Horne-Zeilinger state |2I], which possesses peculiar interest in quantum information science. 

To generate the state, the first step is to produce two pairs of two-dimensional hyper- 
entangled photons, as shown in Figure [TK, via two adjacent thin type-II barium borate 
(BBO) crystals [22] ■ The photon pairs can be generated by spontaneous parametric down- 
conversion (SPDC) [23] and entangled in polarization and frequency DOFs [24j. If pump 
photon downconvert in first (second) crystal, the frequency of the signal photon will be Ui 
{002), and hence the Idler photon takes the opposite. In case a pump pulse coming from left 
(right) threads the first (second) set of crystal, a pair of photons will be produced in the 
spatial modes ai and bi (ci and di). Similarly, a pair of hyperentangled photons can be 
produced in spatial modes 02 and 62 (c2 and ^2) [5J while the pump pulse is reflected by the 
mirror behind the crystal and then threads the crystal second time. Since one can simply 
adjust the relative phase between the first and second possibilities, two pairs of photons in 
the following form are readily obtained: 



Next, we employ the two non-entangled photons a and c in above two pairs to make all 
photons entangled in polarization and spatial mode. Obviously, photons a and c are in the 
state of the form 



{\HH) + \VV)) 0i\cuicu2) + \cu2cu1)) (g)(|ai6i) + 10262)) 



(2) 



{\HH) + \VV)) ^i\uiU2) + \uJ2c01)) (g)(|cidi) + \c2d2)) . 



(3) 



[(1^) + 1^)) (8)(ki) + 1«2))] (S)[m + \v)) (g)(|ci) + |C2))] . 



(4) 
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As shown in Figure [Tt, photons in paths ai and ci (02 and C2) frontally enter Polarizing 
Beam Sphtters (PBSs). Since the PBS transmits H and reflects V polarization photons, we 
can then come to a conclusion that if the two photons pass through paths a[ and c[ and 
C2) simultaneously, it implies that both of the photons are H or V polarized. Afterwards, 
photons in modes and C2 are led to a cross-kerr nonlinear medium [251 126] as shown 

in Figure [21 which brings forth an adjustable phase shift to the coherent state \a) through 
cross-phase modulation (XPM), e.g. 3^, —26, 26 and —6 for four modes, respectively. Then 
we perform the X homodyne detection on the coherent state \a) [27]. In case the measured 
phase shift is 6, the state will be kept. In the end, the reserved states of photons a and c 
can only come out from a[ and c[ or a'g and €3. That means the states of photons a and c 
are in the following entangled form: 

{\HH) + \VV))(^{\a[c[) + \a',c'2)). (5) 

Similarly, we can use another pair of non-entangled photons b and d to make all photons 
entangled in spacial mode and frequency. As mentioned in above, the initial photons b and 
d are in the state: 

[(l^i) + 1^2)) (g)(|&i) + I&2))] (g)[(ki) + \co,)) (g)(Mi) + \d,))] . (6) 

Connecting paths bi, 62? di and d2 with four optical demultiplexers (OD) respectively, which 
are part of Wavelength Division Multiplex (WDM) [28], we then separate each mode into 
two in terms of frequency as shown in Figure [It. The modes bu, 621, du and ^21 correspond 
to frequency coi, while the modes 612, &22, '^12 and ^22 correspond to 1^2 • Taking the same 
procedure as for modes a and c, we lead these eight paths to another cross-kerr nonlinear 
medium, as shown Figure [21 which induces the coherent state \a') phase shifts of 26, 46, —6, 
—36, 36, 56, —26 and —46, respectively. After performing the X homodyne detection on 
\a'), we retain those states with phase shift 6. Then, photons b and d should pass through 
in one pair of the four paths, i.e. bu-du, &i2-'^i2; ^2i-c^2i; and 622-^^22- At last, let those two 
paths originated from the same path via OD link to an optical multiplexer (OM), which is 
other part of WDM and can merge them into one path. For instance, paths &11-&12, which 
both come from bi by OD, are merged into b'l after passing through OM. After taking above 
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procedures, photons b and d then he in the state of the following entangled form: 



iluitui) + \tu2tu2)) 0i\b[d[) + \b',d',)) . 



(7) 



In all, from (EJ, ([3]), ([5]) and ([7]), we find that a four-photon 12-qubit hyperentangled state 
(II]) is created. Unfortunately, the setups shown in Figures [2] and [3] are generally impossible 
to change the sign of the phase shift in single nonlinear medium. The signs of the phase 
shifts must be identical, since they are all proportional to the nonlinear susceptibility [29] . 
To circumvent this difficulty, at least one can use either the cross-kerr with two materials 
of opposite nonlinearities instead of the single material, or use double XPM method [30] 
without minus phase as shown in Figures H] and [51 In Figure HI we use a 50:50 beam splitter 
(BS) to divide the coherent state into two beams \a) \a), and then they are coupled to the 
photonic modes a[ and a'g, c'l and through the XPM respectively. Correspondingly, the 
phase shifts induced by the couplings are 6 and 26 in both beams. Afterwards, the two 
coherent states after the interaction are compared with a 50-50 BS. Project the onto 
the upper beam, in case n = 0, the state is obtained. Analogously, we can obtain the 
state ([7]) with the same procedure as presented in Figure O 

A setup for measuring the hyperentanglement in polarization, frequency and spatial mode 
simultaneously and independently is schematically shown in Figure [Hi Place eight of these 
setups in all eight output paths, the qubits in spatial mode can then be determined. Of 
each spatial mode, the optical demultiplexer splits it into two paths according to frequency, 
which means one therefore obtains the qubits in frequency. Subsequently, at both of these 
two paths the conventional polarization analysis [21] is implemented, by which the qubits in 
polarization are read out. 

3 Quantum characteristics of the hyperentangled state 

To characterize the hyperentangled state ([T]), we evaluate the Bell- type inequality violation 
of it. Note that the 12-qubit state ([T]) is the eigenstate of 2^^ stabilizing operators with unit 
eigenvalue [321 132] • i-e.. 



S,\ij) = \ij),z 



l,2,---2 



12 
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where 

^. = (g)0f (9) 

fc=i 

describe perfect correlations in the state. Here, Oj'^^ll, ±Xj''\ ±Yj''\ -LZ'^p] represent the 
observables for each photon. / and x'^j'\ Z^j^^ denote the identity and Pauh matrices, 
with subscript j = 1,2,3,4 signifying photons a,b,c and d, and superscript k = 1,2,3 
signifying the three DOF. The operators X, Y, and Z satisfy relationship ZX = ~XZ = iV. 

According to Ref. [31] the Bell operator can be expressed as 

2i2 

B = J2S^■ (10) 

i 

The local hidden variable theory(LHVT) sets the upper bound for the expectation value of 
B, and hence satisfies the following inequality [35] 

{B) <maxLHVT\{B)\. (11) 

We notice that among all 4096 Si, there are 1854 of them with negative signs in front as 
shown in the Appendix, like 

In LHVT, the observables are elements of reality, that means their values are precisely 
defined, i.e. either +1 or —1 [36]. Therefore, the upper bound of the expectation value of B 
in (fTOj) is 2^^ — 2 x 1854 = 388. On the other hand, in quantum mechanics, the expectation 
value of B is 2^^ = 4096. Obviously, the QM result greatly violates the inequality f lTTj) . the 
bound set by the LHVT. 

Following, we compare the degrees of violation of different GHZ-type states with same 
amount of qubits but different numbers of photons. According to the Appendix, one may 
conclude that the more the entangled degrees of freedom are, the stronger the violation of 
Bell-type inequality will be. To give a concrete example, in Tabled] we compare the violation 
degrees to the Bell-type inequality of three different 12-qubit entangled states, that is ([T|), 

l^:-') = {\HHHHHH) + \VVVVVV)) (^{\ujiUJ2UJiUJ2UJiUJ2) + \uJ2UJ1UJ2UJ1UJ2UJ1)) , (12) 
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and 

= {\HHHHHHHHHHHH) + \VVVVVVVVVVVV)) . (13) 

The results in the Table indicate that the state ([T]) entangled in three DOFs violates the 
Bell inequality ( ITTl) most, the state (fT2ll entangled in two DOFs goes next, and the state 
(fT3ll entangled in only one DOF violates the Bell inequality least. 

4 Conclusions 

In summary, in this work we propose a novel scheme for the generation of simultaneous 
entanglement of four photons in three independent DOFs, i.e., a 12-qubit hyperentangled 
state, by virtue of the linear optical instrument and cross-Kerr nonlinearity. We show that 
the proposed scheme is within the reach of nowadays technology. Theoretically, our method 
can be simply expanded to the production of hyperentangled state with more photons, 
although practically the difficulties in the manipulation of six or more photons make the 
expansion non-realistic. To check the nonlocality nature of this 12-qubit hyperentangled 
state, we evaluate its violation degree to a Bell-type inequality. The result agrees well with 
the QM prediction while extremely contradicts to the LR constraint, which indicates that 
the genuine multiparticle hyperentangled state possesses certain superiority in comparison 
with the normal entangled state in the verification of QM nonlocality, quantum computation, 
and quantum information processing. 
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Appendix: Comparation of Local Realism Violation 

Here, we compare the violation degrees to the Bell inequality of the GHZ-type states: 
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1^2) = {\HiH2 ■■■//„) + 1^1 ■■ ■ K))® ' , 

\4js) = {\H,H2---Hi) + \V,V2---Vi)), 

where the superscripts 3, 2, 1 represent the number of entangled DOFs; m, n, I signify the 
photon numbers, and satisfy m x 3=^ x 2=1 x 1, i.e they have the same qubits. According 
to the method in calculating the state ([1]) in the text, we just need to calculate the number 
of Si with negative sign, labeled as C{ip). The larger C{ip) means a bigger violation. 

First, we compute the C{ip) of the state which has m photons in one DOF by taking the 
advantage of equation: 

C{ij)=Cl + Cl + ... + C'^'^ (Al) 
with 2i + 4 ^ m. For different values of m, the above equation tells 



( 2*^-2 if m = 4x - 2 

(_1)^2'^ +2'"-2 ifm = 4x-l 

(_l)^2^+2— 2 ifm = 4x 

[{-1)^2'^ +2"^-^ i{m = 4x + l 



where a; = 1, 2 • • ■ . 

Then, we arrive at the results in following ten different cases: 

1. m is odd and m = 4x — 1. In this case n = ^ is not integer, / = 3m = 12x — 3, so 

C(^l)-C(^3) =3((-l)^22--2+24-3)(24-l-(-l)^22-2-24-3)2 + ((_l)^22--2^ 

24x-3)3_((_^)3x26x-3^2i2^-5) ^ _3 x 2-7+6^' X (8 X (- 1)^ + (-l)^' X 2i+^^'+4i+^ -64^) > 
0, for X = 1,2,-- - . 

2. m is odd and m = 4a; + 1. In this case n = ^ is not integer, I = 3m = 12a; + 3, so 
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24a;-i)3 _ (^i^_iyx+226x ^ 212^+1) = -3 X 2-^+^^ X ((-16)^ + (-1)^ + 4^ - 64^) > for 
x = l,2,---. 

3. m is even and m — Sx + 2. In this case n — ^ — 12x + 3, so 

C(V'i) - C(^2) = 3(2S")(28^+2 - 28^)2 + (2^^f - 2((-l)3-+226=^ + 2i2a=+i)(2i2x+3 _ 

(_l)3^+226x _ 2l2x+l^) ^ 2l+12^ + 4^+^2x _ ^ gl+6x > Q for x = 0, 1, • • • . 

4. m is even and m — Sx — 2. In this case n — ^ — 12x — 3, so 

C(^i) - C(V^2) = 3(28^-4) (28^-2 _ 28x-4)2 ^ (28^-4)3 - 2((-l)3^26^-3 + 212^-5) (2^2^-3 - 

(^_;L)3^26^-3 - 2l2x-5^ ^ 4-5+6X ^ (32 + - (_l)^42+3a;) > Q for X = 1, 2, • • • . 

5. m is even and m — 8x + A. In this case n—^ — 12x + 6, so 

C{^l)-C{^2) = 3((-l)2-+224-+l + 2«-+2)(28-+4-(-l)2-+224x-+l_28^+2)2^((_l)2x+224x+l^ 

2^^+2)3 _ 2(2^2^+4) (2120^+6 _ 2i2x+4^ ^ 2^+12^ X (1 - 3 X 2^+^^ + 3 X 4^+^^ + 8^+^^) > 
for X = 0, 1, • • • . 

6. m is even and m = 8x. In this case n = ^ = 12x, so 

C(V'i)-C(V'2) = 3((-l)2^+i24^-i+2«^-2)(28^-(-l)2^+i24^-i-28^-2)2+((-l)2^+i24^-i+ 

28x-2)3 _ 2((_l)3r.+l26x-l ^ 2l2=<^-2) (2l2=^ - (_l)3x+l26r.-l _ 2l2x-2) ^ lQ-l+4x ^ (_;^2 + 

(-1)^ X 2^+2^ - 3 X 2^+4^ + 256^) > for x = 1, 2, • • • . 

7. n is even and n = 4x — 2. In this case Z 2n 8x — 4, so 

C(V'2)-C(V'3) = 2(24^-^)(24^-2-24^-4-((-l)2^24^-3+2«^-6) = 2"^+^^ x (-16+16^) > 
for X = 2, 3, • • • . 

8. n is even and n — 4x. In this case I — 2n — 8x, so 

C(V'2)-C(V'3) = 2((-l)^+i22^-i+24^-2)(24^-(-l)^+i22^-i-24^-2)-((-l)2^+i24^-i + 
28^-2) = 8-^+2^ X (-4 X (-1)^ + 4^) > for X = 1, 2, • • • . 

9. n is odd and n — Ax — 1. In this case Z = 2n = 8x — 2, so 

C(^2) - C(^3) = 2((-l)^+l22^-2 + 24^-3)(24x-l _ (_l)^+l22x-2 _ 24^-3) _ (28^-4) ^ 

2-5+4X X (_4 _ (^_iy X ^ iQx^ > for x = 1, 2, • ■ ■ . 

12 



10. n is odd and n — Ax -\- 1. In this case I — 2n — %x -\-2^ so 

C{i)2) - C(^3) = 2((-l)^+i22^-i + 24^-i)(24^+ii - (-1)^+122^-1 - 24^-1) - (28^) = 
2-i+ix X (_i _ ^_iY X 2^+2^ + 16^) > for X = 1, 2, • • • . 

In conclusion, it is evident that C{ipi) > C{il)2) > Ci'ips), which teUs that for the three 
GHZ-type entangled states with same qubits, the one entangled in three DOFs, the ipi, 
violates LR most, the state ■02 entangled in two DOFs violates the second, and the state ■^s 
entangled only in one DOF violates the least. 
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Table 1: Comparation of violation degree for different entangled states. 



state 


{B)qm 


maxLHVT\{B) 


D = {B)QM/maxLHVT\{B) 




4096 


388 


10.56 




4096 


1024 


4.00 




4096 


1984 


2.06 
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Figure 1: Preparation scheme, a) SPDC source. The two sets of crystal system generate 
photon pairs entangled in polarization, frequency and spatial mode respectively, b) Design of 
constructing the entanglement of four photons with one another in polarization and spatial 
mode. Two PBSs are exploited to implement the selection of polarization, and a cross-kerr 
nonlinear medium which can accomplish the function as quantum nondemolition detector 
(QND) is used to choose photons' spatial modes, c) Design of constructing the entanglement 
of four photons with one another in frequency and spatial mode. Analogously, the selection 
of frequency is achieved by four sets of Wavelength Division Multiplexer (WDM) which 
consists of optical demultiplexer (OD) and optical multiplexer (OM). Moreover, the complete 
experimental set still needs one more cross-kerr nonlinear medium to choose photons' spatial 
modes the same as explained in above. 
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Figure 2: QNDi. The realization of selecting paths of photons a and c with cross-kerr 
nonlinear medium. 
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Figure 3: QND2. The realization of selecting paths of photons b and d with cross-kerr 
nonlinear medium. 
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Figure 4: Double XPM method to realize the function of Fig. |2j 
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Figure 6: Measurement setup. OD: optical demultiplexer; QWP: quarter-wave plate; HWP: 
half-wave plate; PBS: polarizing beam splitter; D: single-photon detector. 
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